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Abstract. A previously-developed colissional, interhemi- 
spheric flux tube model for photoclectrons (PE) [Khazanov et 
al., 1994] has been extended to three dimensions by including 
transport due to ExB and magnetic gradient-curvature drifts. 
Using this model, initial calculations of the high-energy (>50 
eV) PE distribution as a function of time, energy, pitch angle, 
and spatial location in the equatorial plane, are reported for 
conditions of low geomagnetic activity. To explore both the 
dynamic and steady behaviors of the model, the simulation 
starts with the abrupt onset of photoelectron excitation, and is 
followed to steady state conditions. The results illustrate sev- 
eral features of the interaction of photoelectrons with typical 
magnetospheric plasmas and fields, including collisional dif- 
fusion of photoelectons in pitch angle with flux tube filling, 
diurnal intensity and pitch angle asymmetries introduced by 
directional sunlight, and energization of the photoelectron 
distribution in the evening sector. Cross-field drift is shown 
to have a long time scale, taking 12 to 24 hours to reach a 
steady state distribution. Future applications of the model are 
briefly outlined. 


1. Introduction 

An important aspect of photoelectrons (PE) is their trans- 
port through the plasmasphere. The PE escaping from the 
ionosphere experience small-angle scattering when traveling 
through the plasmasphere as a result of the Coulomb interac- 
tion with the thermal plasma. Due to this scattering, some of 
the PE are scattered to the outside of the loss cone and undergo 
magnetic reflection, i.e. they become trapped (for a compre- 
hensive review of this problem see Khazanov et. ai [1993- 
1995]). 

Among the numerous papers devoted to PE transport be- 
tween the magnetoconjugate regions of the ionosphere, no re- 
sults have yet been published which describe the evolution and 
formation of PE fluxes under nonstationary conditions on a 
global scale. The previous models just picked a single field 
tube, or part of one, and considered the PE transport problem 
only inside this restricted region. However, PE calculations 
are an important part of global modeling of the plasma distri- 
bution in the inner magnetosphere. The authors believe that 
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the investigations of non-steady-state PE transport on a 
global scale is the logical next step on the path to developing 
a self-consistent model of the thermal plasma in the inner 
magnetosphere and can form the basis of a numerical space 
weather prediction scheme. Such a study also has particular 
relevance to the goals of the IMI (Inner Magnetosphere 
Imaging) mission currently being developed at the Marshall 
Space Flight Center. The temporal evolution and spatial dis- 
tribution of plasma parameters in the inner magnetosphere and 
the variation in these quantities with changing solar and mag- 
netic activity arc critical factors in setting the sensitivity 
ranges and spatial and temporal resolution requirements for 
IMJ/ENA (Energetic Neutral Atom) and EUV (Extreme 
Ultraviolet) imaging instrumentation. 

In this paper, the first results of a numerical study of the 
non-steady-state bounce-averaged kinetic equation for PE, in- 
cluding their transport due to ExB and magnetic gradient- 
curvature drifts in the inner magnetosphere, are presented. The 
distribution function in time, spatial location in the equatorial 
plane, energy, and pitch angle are among the parameters calcu- 
lated by the model. To insure the applicability of the bounce- 
averaged equations, the present study considers electrons with 
energies greater than 50 eV. 


2. Model 


The two main processes controlling the behavior of PE in 
the plasmasphere are their motion in the inhomogeneous 
geomagnetic field and Coulomb collisions with the back- 
ground thermal plasma. The kinetic equation for PE in the 
guiding center approximation at altitudes greater than about 
800 km can be written as [Khazanov et ai , 1994]: 
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where {t, s, R lt p, v, E} are time, distance along a given 
field line, position vector in the equatorial plane, the cosine 
of the pitch angle, velocity, and energy, respectively; 
v il=pv is the velocity along the magnetic field; V D is the 
drift velocity due to ExB and magnetic gradient-curvature 
drifts; E and B are the electric and magnetic field intensities, 
respectively; / is the distribution function of the electrons; 
A=27te 4 lnA= 2.6x1 O' 12 cm 2 eV 2 (e is the electron charge, InA 
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is the Coulomb logarithm); and n e is the density of the (her- 
mal plasma. 

Performing the standard procedure of bounce-averaging (1) 
[see, for example, Lyons and Williams , 1984 or Gefan and 
Khazanov. 1990], the number of independent variables can be 
reduced to five and the kinetic equation (1) can be presented in 
the bounce-averaged form: 
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The brackets (<!;) denote bounce-averaging: 
ds 
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where £ = /, V D , n e 


(2) 


and the following notation was introduced: 
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with and denoting the magnetic field and the cosine ot 
the pitch angle at the magnetic equator of the flux tube. 

Note that equation (2) is not applicable to describe the PE 
distribution function in the loss cone. It can only be used for 
the description of trapped particles. The loss cone is defined 
by Bob and lhc trap ping regio n by 

<|ju„|<ju„ fc , where »„ h = Jl-B 0 /B{s,) is the 
loss cone boundary and B(s,) is the magnetic field at the low 
altitude boundary (the region over which the electron distribu- 
tion function is defined in terms ot s and Bo * s shown in Figure 
1 of Kkazanov et al. [1992]). That means, in order to solve 
equation (2), it must be coupled with a corresponding loss 
cone equation [Khazanov et al., 1992] or a boundary condition 
for the electron distribution function at the edge of the loss 
cone should be introduced. 

Equation (2) has been transformed into the conservative 
form of lhc variables (l, R,„<p,E,fX„), where R„ is the radial 
distance in the equatorial plane, and (p is the geomagnetic east 
longitude, and then solved numerically in light of the ring 
current model of Jordanova et al. [1995], which was adjusted 
for the PE energy range and coupled with a field-aligned model 
[Khazanov and Liemohn , 1995; Licmohn and Khazanov , 
1995]. Utilizing the approach developed by Fok et al. [1993], 
equation (2) is decomposed using a double splitting method, 
and advance the solution in a single dimension at each half 
time step. The order of the solution operators is reversed in 
the next half time step to achieve second order accuracy in 
time. The drift terms and the Coulomb drag energy term are 
discretized by the Lax-Wendroff scheme and the first order up- 
wind scheme. The superbee flux-limiter method is used to 
switch between the high-order Lax-Wendroff scheme, which 
works well in smooth regions, and the low-order upwind 
scheme, which behaves well in the presence of strong gradi- 
ents. This high resolution method gives second order accuracy 
on smooth solutions and well-resolved nonoscillatory discon- 
tinuities [Leveque, 1992]. The Crank-Nicholson scheme is 
applied for discretization of the Coulomb pitch angle diffusion 
term according to Khazanov [1979]. 
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Figure 1. L=5 daysidc pitch angle distributions. 


The following boundary conditions were used: a) outside the 
region under consideration, L=1.75 to L=6.5, PE fluxes equal 
zero; b) the energy range is between 50 eV and 200 eV with 
zero flux in the tail of the energy distribution; c) the solution 
is assumed to be periodic in azimuth ( (p ); d) the only source of 
superthermal electrons is photoionization and secondary ion- 
ization due to these particles; and, e) the solution is assumed 
to be symmetric at 90° pitch angle, and during the day the PE 
distribution function at the loss cone boundary has been calcu- 
lated by using the Khazanov and Liemohn [1995] model. 
During the night, the absorption of PE in the ionosphere, due 
to the emptying of the loss cone (twice per bounce period, 
I#), is taken into account by extending (2) into the loss cone 
and introducing here the loss term according to Lyons and 
Williams [1984]: 



A 


where 


{ t b /2 in the loss cone 
oo out of the loss cone 


In this case it is required to use the boundary condition at zero 
pitch angle which is the same as in the Khazanov at al. [1994] 
model: d(f)/dBo = 0. 

When considering PE formation in the plasmasphere on a 
global scale, it is necessary to keep in mind the possibility of 
PE interaction with the vast collection of plasma waves sup- 
ported by the magnetospherie plasma. The wave-particle in- 
teraction mechanisms influencing the PE distribution have not 
been incorporated into this model yet, and as a first step will 
not be included in the calculation. The justification for omit- 
ting this process is as follows. The ring current and radiation 
belts are two major mesoscale phenomena that interact with 
the thermal plasma to generate most of the plasma waves in 
the region under consideration. Ring current-plasmasphere in- 
teraction amplifies Alfven, fast-magnetosonic, and ion-cyclo- 
tron waves, while radiation belt-plasmasphere interaction am- 
plifies whistler mode waves. The former waves resonate with 
electrons at energies up to a few electron volts, and the latter 
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Figure 2. L=5 nightside pitch angle distributions. 


resonate at keV energies. Among these waves, the most ap- 
propriate candidate for additional photoelectron scattering is 
interaction with whistlers through Cherenkov resonance, but 
this will be much less efficient than Coulomb scattering for 
low levels of geomagnetic activity [Kuzivanov and Khazanov, 
1984]. Although the energy range under consideration is usu- 
ally quite stable, it degrades into the low-energy (E<50 eV) 
distribution, which can be unstable and generate plasma waves 
itself Khazanov [1979]. 

The following input was used to perform the PE model calcu- 
lations. The bounce-averaged drift velocities of the guiding 
center can be obtained for the ExB and magnetic gradient- 
curvature drifts, considering that only the equatorial electric 
field is contributing to the bounce-averaged drift [ Roederer , 
1970]. The geomagnetic field in this simulation is taken to be 
a 3-dimensional dipole tield. The electric field is expressed as 
the gradient of the Volland-Stcrn semiempirical potential 
model [Volland, 1973; Stern , 1975]. All of the calculations 
were performed for a quiet geomagnetic condition K p =l (when 
the main superthermal electron source is photoionization) on 
March 21, 1986, with F10.7 and <F10.7> values of 150, 
chosen so the atmospheric conditions are symmetric and the 
solar radiation is at an average intensity level. The plasma- 
spheric cold electrons are calculated with the model of 
Rasmussen et al. [1993]. 

Here, the first results from this new global PE model are 
presented and the processes of ionosphere-plasmasphere- 
magnetosphere PE coupling and distribution function 
formation on a global scale are discussed. 


3. Results and Discussion 

To describe PE formation on a global scale, the global PE 
transport is calculated on the time scale up to 48 hours. The 
early stages of the refilling process and the behavior of super- 
thermal electrons in the geomagnetic trap on the time scale up 
to one hour were considered using a recently developed model 
[Khazanov et. al . , 1993-1995; Liemohn and Khazanov , 


1995], which represents a unified approach by self-consis- 
tently coupling the interaction of the two hemispheres and the 
capture of PE in the plasmasphere. In these papers, the readers 
also can find the typical PE energy and pitch angle distribu- 
tion in the region of the source (below 800 km) that has been 
used as the dayside loss cone boundary in this study. 
Considering PE transport based on the bounce-averaged 
kinetic equation (2), it is necessary to keep in mind that this 
approach is valid since the bounce period of PE is much 
shorter than the scattering processes' decay lifetimes; 
presently just Coulomb collisions. For the time period com- 
parable with or less than the bounce period, the kinetic equa- 
tion for PE in the guiding center approximation (1) should be 
used, using the approach proposed by Khazanov et. al. [1993- 
1995]. 

In this paper, PE results above a low energy limit of 50 eV 
are presented, where the applicability of the bounce-averaged 
kinetic equation (2) in the region under consideration (L=L75 
to L=6.5) is obvious. An extension of (2) to the low energy 
range, where the bounce period of PE is comparable with the 
time scales of their scattering processes, will be investigated. 
This would allow the inclusion of practically all important PE 
energies contributing to magnetospheric thermal plasma en- 
ergy deposition. 

It is interesting to start the calculations with initially 
empty flux tubes at t=0 on the day and night sides of magneto- 
sphere and to consider the non-steady-state development of 
the photoelectron distribution function in velocity and con- 
figuration spaces as a result of photoelectron energy pumping 
and pitch angle redistribution in the geomagnetic trap on a 
global scale. The calculations in this paper have been per- 
formed only for the case of low geomagnetic activity, K p =l; 
however, the model can handle any conditions of geomagnetic 
activity or K p history with corresponding variations in the 
neutral atmosphere, ionospheric and plasmaspheric densities, 
and electric field of magnetospheric convection. 

Figure 1 shows the pitch angle distribution formation on 
the dayside at L=5 for the energies 50 eV, 100 eV, and 150 eV 
in the different magnetic local time (MLT) sectors: MLT=6 is 
morning-sunrise (a-c), MLT=12 is noon (d-f), and MLT=18 is 
evening-sunset (g-i). The different lines indicate the time of 
the solution, which starts with the initial PE distribution 
located only in the loss-cone at t=0 and evolves to t=48 hours. 
The results of the model are shown here in terms of PE fluxes, 
connected with the electron distribution function according to 
<f> = 2Ef /m~ , where m is the electron mass. As seen in the 
presented results, in the noon-evening sector the pitch angle 
distribution has similar pitch angle and energy dependencies 
to those calculated in previous papers [Khazanov et al., 1993- 
1995]. For example, the shape of the pitch angle distribu- 
tions in this time sector remains the same, it also has the same 
energy dependence in the transition region from the small 
pitch angles near the loss cone to the large pitch angles of the 
trap zone, and the time scale of the refilling process is propor- 
tional to E" 1/2 . In the morning, the evolution of the pitch 
angle distribution has a more complicated character and is sim- 
ilar to the dayside pitch angle distribution only during the first 
few hours of plasmaspheric refilling. After 10-12 hours, the 
shape of the pitch angle distribution in the trap zone on the 
morning side changes. The tail of the pitch angle distribution 
is trying to stretch out. There appears to be another source 
besides the loss cone, and the part of the pitch angle distribu- 
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tion that rises depends on the PE energy. To explain such 
unusual pitch angle distribution behavior, attention should be 
given to what is going on in the nightside magnetosphere. 

Figure 2 demonstrates the pitch angle distribution forma- 
tion on the nightside at L=5 for the energies 50 cV, 100 eV, 
and 150 eV in different MLT sectors: MLT=21 is the late 
evening (a-c), MLT=0 is local midnight (d-f), and MLT=3 (g-i) 
is deep night. The different lines indicate the time of the solu- 
tion, which starts with the zero PE distribution function on the 
nightside, and evolves to t=48 hours. The source of the PE on 
the magnetospheric nightside is the drift motion of the par- 
ticles from the dayside due to the ExB drift in the presence of 
the corotation and magnetospheric convection electric fields. 
The main feature of the nightside PE distribution function is 
the sharp depletion of intensity in the loss cone, because after 
MLT=18, photoionization in the ionosphere disappears. PE 
flux in the loss cone does not go to zero, however. Depending 
on MLT and equatorial distance, two processes, Coulomb col- 
lisions and magnetospheric convection, move particles from 
the trap zone to the loss cone. These plots clearly show the 
time-dependent dynamics of the pitch angle development in 
the different MLT sectors. 

As seen in Figure 2 at the deep night sector, the pitch angle 
distribution is flattening out at large pitch angles. This fea- 
ture can be explained as a result of the conservation of the first 
and second adiabatic invariants as the PE drift from the dusk 
sector to the dawn sector along convective trajectories 
moving toward the Earth. Going back to the morning sector 
in Figure 1 (a-c), such unusual pitch angle distribution beha- 
vior in this region can now be explained simply by delivering 
particles from the nightside of the magnetosphere to the 
morning sector. Similar ideas can be proposed to explain the 
difference in pitch angle distribution behavior at MLT=21 (a- 
c), MLT=0 (d-f), and MLT=3 (g-i). In this case, the late 
evening pitch angle distribution retains some of the shape of 
the dayside PE distribution. 

4. Summary 

A global PE model including transport across geomagnetic 
field lines due to ExB and magnetic gradient-curvature drifts 
has been developed. Using this model, a detailed calculation 
was performed of the PE distribution as a function of time, en- 
ergy, pitch angle, and spatial location in the equatorial plane 
on a global scale and presented the first results of this study 
for a low' level of geomagnetic activity. 

The authors believe that the investigation of non-stcady- 
state PE transport on a global scale is the logical next step on 
the path to developing a self-consistent model of the thermal 
plasma in the inner magnetosphere and can form the basis of a 
numerical space weather prediction scheme. Here, only results 
with photoionization and secondary production of superthcr- 
mal electrons are presented, but the model is capable of hand- 
ling various other sources, such as energetic precipitation and 
sunward How from the magnetotail. The ultimate objective of 
this global superthermal electron model will be to couple it 
with the thermal plasma and the ring current and predict the 
energy deposition from superthermal and hot populations to 
the cold magnetospheric plasma. Extension of the energy 
range down to a few eV and up to several tens of keV is also 
planned and will provide the numerical forecast of the energy 


deposition from the superthermal electrons to the thermal 
plasma and the hot electron population on a global scale. 
This would allow the use of these results in cold magneto- 
spheric plasma models to predict plasma temperatures and den- 
sities of all significant thermal plasma species on a global 
scale. 
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